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In 1], the temperature dependence of the probability density function (PDF) of the magnetization (the order 
parameter) was first shown by numerical Monte Carlo analysis using the finite 2D XY model. This result contradicted 
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. the previous claims found in the literature 0. Recently, shows semi-analytically that the PDF depends on the 
^\ ' system temperature T . To show this dependence the authors compute the contribution of the multi loop graphs 
(MLG) to the moments of the PDF within the Harmonic approximation of the XY model. They show that the MLG 
depend on T and may not be neglected in the thermodynamic limit contrary to previous assumptions 2j. Using a 
Monte Carlo simulation they compute the skewness 73 or normalized third moment of the PDF, which "provides 
a clear measure of the variation of the PFD with temperature", and using a Monte Carlo simulation, obtain the 
numerically approximated expression for 73, valid for a square lattice of lattice size L = 16, 
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73 (T) w -0.85 + 0.126T - 0.0048T 
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They also obtained for the lattice size L = 32 "a value much closer to the theory" with 73 (T) w —0.88 + 0.15T, 
and finally argued that the skewness is relatively computational expensive due to the need for averaging and, as their 
results appear to confirm the evolution of 73 with T, they "... leave the determination of the precise form of 73 (T) 
I from larger systems to another time" . 

This comment is devoted first to show that an explicit analytical expression for 73 (T), valid for arbitrary system 
size and to all orders in T, can be deduced starting from their equation (11) for the higher moments (m p ) of the 
PDF, and second to show that the numerical values of the slope of 73 reported in their paper are not just numerically 
incorrect, but their scaling with the lattice size is wrong. 
fN| \ Starting from their equation (11) and using the translational invariance of the lattice propagator G, the exact 
^ > expression for the moments to all orders in T can be obtained, yielding: 
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VO- (M)=exp[-TG(0)/2] (1) 
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( M 2 ) = ]T cosh[T(G(l?)] (2) 

{ M 3 ) = ^ E { exp(-TG(^)) cosh[T(G(^) + G(l? - ?)] + exp(TG(^)) cosh[r(G(l/) - Gtf - V)] J (3) 



where A denotes the lattice, N — L 2 is the volume and the lattice propagator G is given for example by its Fourier 
representation 
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Kl is the lattice momentum defined as usual as (ifi)i = 2sin(.Ki/2) with i = 1,2 and Ki lies in the first Brillouin 
zone. We now use the definition of the skewness as the third normalized moment 73 (T) = / \(M — ( M ) ) / a ] 3 ) 
and write 

73 (T) = [( M 3 - 3 ( M 2 ) (M} + 2( M ) 3 )] (4) 
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and explicitly expanding up to order T, we obtain 

73(T)^. 93 (l) 3/2 {l-^T + (T 2 )} (5) 

where the quantities g n are defined in terms of the power n of the lattice propagator G as g n = G n (0)/N n ~ 1 . The 
expression of eqn. JSJ) agrees with the corresponding equation (33) of ref. Q, where the skewness and kurtosis are 
computed for the full 2D XY-model up to two-loops, including the anharmonic corrections to the Hamiltonian (see the 
two last terms in eqn. (33)), which are suppressed by a volume factor N, and which are negligible in the thermodynamic 
limit but are relevant for finite lattice sizes. Higher order corrections in T can be directly computed along the lines 
outlined here. From the numerical point of view, we use the known numerical values for the lattice coefficients g n , 
which are directly and easily evaluated by using MATLAB for example, and obtain 73 (T) « —0.8540 + 0.1358T for 
L = 16 and 73(7) ~ —0.8763 + 0.1331T for L = 32 respectively. Independent of the numerical differences found for 
the first two coefficients of the skewness, one observes that the values for the slope (linear term in T) reported in Q 
increases with the system size, contrary to the tendency of the analytic expression displayed by eqn. JSJ. In fact, the 
slope is a decreasing function of the system size and its thermodynamic limit converges to the value 0.1319. Finally, 
and in order to compare our results with the accurate numerical values obtained in 5] for the skewness, we evaluate 
the contribution obtained from the anharmonic corrections 
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and obtain Sj 3 (T) w 0.0197T for L = 16 and Sj 3 (T) w 0.0191T for L = 32 respectively . Perfect agreement is found 
with the corresponding value in (see eqn.( 36)), when adding this contribution to the value obtained in eqn. © 
for 73 (T). 
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